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We investigate the future evolution of the dark energy universe in modified gravities including 
F[R) gravity, string-inspired scalar-Gauss-Bonnet and modified Gauss-Bonnet ones, and ideal fluid 
njrs', with the inhomogeneous equation of state (EoS). Modified Friedmann- Robertson- Walker (FRW) 

^— s . dynamics for all these theories may be presented in universal form by using the effective ideal fluid 

("■^ ' with an inhomogeneous EoS without specifying its explicit form. We construct several examples 

^sj , of the modified gravity which produces accelerating cosmologies ending at the finite-time future 

singularity of all four known types by applying the reconstruction program. Some scenarios to 
^ ' resolve the finite-time future singularity are presented. Among these scenarios, the most natural one 

^S \ is related with additional modification of the gravitational action in the early universe. In addition, 

late-time cosmology in the non-minimal Maxwell-Einstein theory is considered. We investigate the 
forms of the non-minimal gravitational coupling which generates the finite-time future singularities 
and the general conditions for this coupling in order that the finite-time future singularities cannot 
emerge. Furthermore, it is shown that the non-minimal gravitational coupling can remove the 
finite-time future singularities or make the singularity stronger (or weaker) in modifled gravity. 
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INTRODUCTION 



^ . Recent observational data strongly indicate the existence of the dark energy, which generates the accelerating 

l/~) ' expansion of the present universe. In particular, the five-year Wilkinson Microwave Anisotropy Probe (WMAP) 
t — ■ ■ data [11] give the bounds to the value of the equation of state (EoS) parameter wde, which is the ratio of the pressure 
*/~J ■ of the dark energy to the energy density of it, in the range of —1.11 < iude < —0.86. This could be consistent if the 
^^ , dark energy is a cosmological constant with ujde = ^1 and therefore our universe seems to approach to asymptotically 
[ — ' de Sitter universe. It is also believed that there existed the period of another accelerating expansion of the universe, 
^^ I called inflation, in the early universe. In many models of inflation, the accelerating expansion could be generated 
00 . by almost flat potential of scalar field(s), called inflaton. Hence, in the period of the inflation the universe could be 
described by (almost) de Sitter space. Thus, there is the striking similarity between the very early and very late 
universe. 

Although both of the accelerating expansions seem to be de Sitter type, the possibility that the current acceleration 
^ , could be quintessence type, in which wde > —1, or phantom type, in which ujde < —1, is not completely excluded, 
j^ ■ Furthermore, even if the current accelerating universe is described by ACDM epoch, it is quite possible that it may 
enter to quintessence/phantom phase in future. Similarly, even if the early-time inflation may be de Sitter type, there 
may exist pre- inflationary stage in which the evolution of the universe is different from de Sitter type, e.g., it could be 
quintessence/phantom epoch. It is often assumed that the early universe started from the singular point often called 
Big Bang. However, if the current (or future) universe enters the quintessence/phantom stage, it may evolve to the 
finite-time future singularity depending on the specific model under consideration and the value of the effective EoS 
parameter. This suggests that the same theory should describe the whole history of the expansion of the universe. 
To achieve this, two scalar fields, one field (inflaton) to describe the inflation and another field to do the dark energy 
have often been introduced. 

In this paper, we consider another approach, which might be more natural than introducing two scalar fields in order 
to unify the early-time inflation with late-time acceleration: that is modified gravity (for a review, see [3])- In this 
approach, one starts from some unknown fundamental gravity. At the very early universe in which the curvature is very 
large but quantum gravity effects may be neglected, the restricted specific sector of such a theory predicts inflation. 



Also at Lab. Fundam. Study, Tomsk State Pedagogical University, Tomsk 



In course of the evolution, the curvature decreases and the next-to-leading terms become relevant at the intermediate 
universe (the radiation/matter-dominated stage). Note that the observational data seems to be consistent if the 
intermediate universe could have been governed by the standard general relativity. As the curvature becomes smaller, 
the universe enters the dark energy epoch controlled by the different sector of the unknown fundamental gravitational 
theory different from general relativity. The corresponding gravitational terms are leading ones in comparison with 
the ones of general relativity at the current curvature. Hence, the evolution of the universe defines the modified 
gravitational theory predicting its evolution at each stage [3|. From other side, the effective evolution of modified 
gravity is responsible for the history of the expansion of the universe. This also indicates that a right approach to 
the understanding of the fundamental gravitational theory is to study the history of the expansion of the universe, 
which will give the information about the leading sectors of modified gravity at each epoch. Moreover, the consistent 
examples of modified gravitational theories passingthe local tests and unifying the early-time inflation with late- 
time acceleration have already been constructed [J, l5| . To understand such theories more clearly, it is reasonable to 
investigate these theories at extreme situations, for instance, near to singular points where some fundamental features 
of the theories may be discovered. 

Besides dark energy, there is another unknown component in the universe, namely, dark matter. While dark energy 
has the large negative pressure, the pressure of dark matter is negligible. In many scenarios, dark matter could originate 
from particle physics, e.g., dark matter could be a lightest supersymmetric (SUSY) particle, or particle coming from 
extra dimensional theories. There exists another kind of scenario in which dark matter would be explained by the 
modification of the gravity, as in case of dark energy considered in this paper. There are many such scenarios like a 
modified Newtonian dynamics (MOND) Q, Tensor- Vector-Scalar gravity (TeVeS) Q, and the generalized Einstein- 
Aether theory Q derived from the Einstein- Aether theory [0| and the equivalent one in non-relativistic limit to MOND. 
Dark matter might be explained by F{R) gravity [i3|, where F{R) is an arbitrary function of the scalar curvature R. 
There is a scenario in which dark matter particles could be given by the scalar mode in F{R) gravity (the first paper 
in [Hi). 

In the present paper, we study the future evolution of the dark energy epoch for various modified gravities: F{R) 
gravity, scalar-Gauss-Bonnet and modified Gauss-Bonnet ones, and effective ideal fluid with the inhomogeneous EoS, 
which includes the explicit dependence from Hubble rate and curvature. In particular, we are interested in the behavior 
of the accelerating cosmological solutions in these theories when those solutions approach to the finite-time future 
singularity. Of course, not all of these theories predict such singularities. It depends on the value of the effective EoS 
parameter and the structure of theories. For instance, the expansion in the phantom phase which is not transient 
predicts a future Big Rip singularity. 

The paper is organized as follows. In Sec. II, we present the modified Friedmann-Robertson- Walker (FRW) dynamics 
in a universal way. The following models are considered: F{R) gravity, scalar-Gauss-Bonnet and modified Gauss- 
Bonnet theories as well as ideal fluid with the inhomogeneous EoS. It is indicated briefly how the history of the 
expansion of the universe can be reconstructed for such a universal formulation. Sec. lIIIl is devoted to the study of the 
finite-time future singularities in F(R) gravity. Using the reconstruction technique, we present several examples which 
predict the accelerating FRW solutions ending at the finite-time future singularity. It is demonstrated that not only the 
Big Rip but other three types of the finite-time future singularities may appear. In Sec. lIVl we discuss various scenarios 
to resolve the finite-time future singularities. The most natural scenario is based on the additional modification of 
the inhomogeneous EoS or the gravitational action by the term which is not relevant currently. However, such a 
term which may be relevant at very early or very late universe may resolve the finite-time future singularities. It 
is interesting that the presence of such a next-to-leading order term does not conflict with the known local tests. 
Another scenario is related to the account of quantum effects which become relevant near to the flnite-time future 
singularities. Sec.|V]is devoted to the construction of scalar-Gauss-Bonnet and modified Gauss-Bonnet gravities which 
predict the late-time acceleration ending in the finite-time future singularities. Using the reconstruction technique, 
we present the corresponding effective potentials. 

In Sees. IVll to IVIIll we study some cosmological effects in the non-minimal Maxwell-Einstein gravity with general 
gravitational coupling. In Sec. IVIl we describe our model and derive the effective energy density and pressure of the 
universe. In Sec. IVIIl we consider the finite-time future singularities in non-minimal Maxwell-Einstein gravity. We 
investigate the forms of the non- minimal gravitational coupling of the electromagnetic field generating the finite-time 
future singularities and the general conditions for the non-minimal gravitational coupling of the electromagnetic field 
in order that the finite-time future singularities cannot emerge. Furthermore, in Sec. IVIIll we consider the infiuence 
of non-minimal gravitational coupling of the electromagnetic field on the finite-time future singularities in modified 
gravity. It is shown that a non-minimal gravitational coupling of the electromagnetic field can remove the finite-time 
future singularities or make the singularity stronger (or weaker). Some summary and outlook are given in Sec. IIXI 
We use units in which k^ = c = h = 1 and denote the gravitational constant SttG by k^, so that k^ = Stt/Mpi^, where 
Mpi = C?~^/^ = 1.2 X lO^^GeV is the Planck mass. Moreover, in terms of electromagnetism we adopt Heaviside-Lorentz 
units. 



II. MODIFIED FRW DYNAMICS AND THE RECONSTRUCTION OF THE HISTORY OF THE 

EXPANSION OF THE UNIVERSE 

In this section, we present the general point of view to the modification of the FRW equations which may be caused 
by alternative gravity or ideal fluid with the complicated EoS. We consider F{R) gravity, scalar-Gauss-Bonnet and 
modified Gauss-Bonnet theories, and ideal fluid with the inhomogeneous EoS. In addition, we briefly indicate how 
the history of the expansion of the universe can be reconstructed through such a universal formulation. 

The flat FRW space-time is described by the metric 

ds^ ^-dt^ +a^{t)dx^ , (2.1) 

where a{t) is the scale factor. In the Einstein gravity, the FRW equations are given by 

P^^H\ p^-^{2H + iH^) , (2.2) 

where H = a/a is the Hubble parameter, a dot denotes a time derivative, ' = d/dt, and p and p are the energy 
density and pressure of the universe, respectively. We have assumed the flat three-dimensional metric in accord 
with observational data. Let us consider any modified gravity (for a review, see [2]) like F{R) gravity (for reviews, 
see [3, [ill), the scalar-Gauss-Bonnet one, or the modified Gauss-Bonnet one {F{Q) gravity, where Q is the Gauss- 
Bonnet invariant given by Eq. (|2.f6p below). In this case, the part of modified gravity may be formally included 
into the total effective energy density and the pressure as in Ref. |12l |. in which general inhomogeneous EoS fluid is 
introduced. In this case, the modified FRW equations have the well-known form: 

Peff = 4^' ' Pcff = -^ hil + 3hA , (2.3) 

where pcff and pcs are the effective energy density and pressure of the universe, respectively. Note that the contribution 
of (unusual) ideal fluid should also be included in the left-hand side (l.h.s.) of the above modified FRW equations p^ . 
Poff and Pcff satisfy more general EoS as 

Pcff = -Pcff + / (Pcff) + G (i/, H, if , • • • ) , (2.4) 

or even more complicated one. The other point of view is possible: one can only keep the contribution of matter in 
the energy density, while the gravitational modification should be parameterized by the above function G. We should 
note that peff and pcff defined in (j2.3p satisfy the conservation law identically: 

Poff + 3i7(peff+Peff)=0 . (2.5) 

As an example, we may consider the F{R) gravity [3, lll| whose action is given by 



Sfir) = / d^x^^ -V + ^m , (2.6) 



where F{R) is a proper function of the scalar curvature R and Cm is the matter Lagrangian. One may separate the 
modified part in F{R) from the Einstcin-Hilbert one as 

F{R)^R + f{R). (2.7) 

In the FRW background with flat spatial part, pcs and Pcs are given by 

Peff = -^ {'l-f^^^ + 3 (if2 + h) f'(^R) _ 18 (^iH^H + HH) f"{R)\ + p„,attcr , (2.8) 

Pcff = ^ ( ^/(^) - (3-^' + ^) /'(-R) + 6 (Si/'ii + 4ii2 + QHH + 'ii) .f"{R) + 36 (aHH + H^ ' J"'{R) 

+Pmatter ■ (2.9) 

Here, Pmatter and Pmatter are the energy density and pressure of the matter, respectively, and the scalar curvature R 
is given by i? = 12H^ + 6H. If the matter has a constant EoS parameter w, Eq. (|2.4p has the following form: 



Pcff = wpcff + G {H,H,H,- 
I /I 



g[h,H,---) = -(-J^f{R)-{3{l+w)H^ + (l + 3w)H},nR) 

+qU% + 12w)H'^H + AH'^ + {& + ?.w)HH + h\ f'\R)+mUHH + H\ f"{RU . (2.10) 



Let us consider several examples. In case of the model [1^ in which f(R) is given by 



f{R) = -j^+PR" 



(2.11) 



one finds 



G (h, H, 



1 / 1 + w / a 



^ + PR"') -{3(1 + w) H^ + {1 + 3w) ij} (-^ + 7i/3i?"-i 



+6 {(8 + 12w) H^H + 4iJ2 + (6 + 3u>) HH + H^( -—^ + n{n ~ 1)I3R 



+36(AHH + H) ( — j+n(7i-l)(n-2)i?"-3 



For the Hu-Sawicki model UA 



fHs{R) 



m?ci [R/m^) rn?ci m^Cxjci 



C2{R/m^f + l C2 C2{R/m?f + l 



(2.12) 



(2.13) 



we get 

G[H,H, 



l + w ( rri^ci [R/m?Y 



2 \c2 {R/m? 
|(8 + 12w) H'^H + 4iJ2 + (6 + 3w) FF + h\ 

-36 UhH + Hj ' "^ ^ ' 



- j + {3 (1 + w) i72 + (1 + -iw) i^} [ - 



nci {Rjm} 



Xc2{Rlm^T^l) 



(C2 {Rlm-^f + 1)' (C2 (i?/TO2)" + 1)3 



2 / "("-^)(r'^'=^ {Rim') "-' ^"'^"-^^^^^^ (i?/m2) '"-' 5!^ (i?/TO2) 3"-3 \ n 



(C2 (i?/m2)" + 1)' 



(C2 (i?/m2)" + 1)' 



(C2 (i?/TO2)" + 1)^ 



(2.14) 



Note that recently the observational bounds for F(R) theories have been discussed in Ref. [Ig- We may also consider 
the F(Ci') gravity [l^, whose action is given by 



5 = 1 d^x./—g I ^ (i? + !q m + C,n I 



where Q is the Gauss-Bonnet invariant: 

g = R'- 4i?^,i?^'^ + i^^.^.i?'"-?" 

In the model, the effective energy density and pressure are given by 



(2.15) 



(2.16) 



PeS 



PoS 



1 

2^ 
1 

2^ 



Pmatter ; 



GfgiG) - fg{Q) - 24' H^ [2H' + HH + AH'Hj f^ 
fg{g) + 24' H' (SH'^ + 20H'H' + QH^ + 4H'^H + H'h) fg{g) 
(2H' + HH + 4H'h) f^'{g) 



" PniSLttcr 



(2.17) 



In the FRW background, we find g ^ 24 ( H'H + if* ) . If we assume that the matter has a constant EoS parameter 
w, again, Eq. (|2.4p has the following form: 



Gg \H,H,- 



Poff = wpos + Gg [H, H, H, 

:^[ii + w)fgig)-wgf^ig) 

r2iV2 , ciV3 



24^ii^ 3iJ** + 2QH'H' + 6H'^ + 4H'^H + H'h + wH^ [2H' + HH + 4H'H) f^{g) 



-24^H^ (2H' + HH + 4H'h] f^'{g) 



(2.18) 



An example is given by [if 



fgiG)^fo\G 



1/2 



(2.19) 



where /o is a constant. In this case, one gets 



Gg [H,H,- ■ ■ ] — 7^-2 



A 

2k2 



w 



1+f iiar^' 



lUH^ ( 3ff^ + 20H^H^ + 6H' + 4H^H + H^H + wH^ [ 2H^ + HH + AH'H 

2 g 



1 



-9 • 2A^H^ Oih'^ + HH + 4H^H 



2\g\ 



7/2 



|g|3/2 

(2.20) 



In the same way, one can obtain the function of modified gravity G for other models including non-local gravity [l7j . 
It is interesting to note that the perturbations of the above theory should be considered by using the analogy with 
effective field theory |18j . 

We now study the general reconstruction program in terms of G in (j2.4p . For simplicity, we assume that the matter 
has a constant EoS parameter w. Using (|2.3p . we find 



G[H,H, 



2H + 3(1 + w)H' 



(2.21) 



Let the cosmology be given by iJ = H{t). The right-hand side (r.h.s.) of ()2.2ip is given by a function of t. If the 
combination oi H, H, H ■ ■ ■ in G (H,H,- ■ ■) reproduces such a function, this kind of cosmology could be realized. 
As an illustrative example, we may consider the case in which H{t) is given by 



H = h, + ^ 



which gives 



H 






H 



2hi 



and the r.h.s. in (|2.2ip is given by 

^ ''2H + 3{l + w)H' 



^ f-if^ 1 ^;2 , Hi + w)hohi 
-2 (3(1 + w)h„ + 



-2hi + 3{1 + w)hl 
f2 



For (trivial) example, if G is given by 

G (h, i^) = 4 {-3(1 + w)hl + 6(1 + w)hoH + [2 - 3 (1 -f «;) h^] ij} , 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



(|2.22p is a solution. Of course, there is large freedom of the choice oi G lH,H, H,- ■ ■ ) , but the form could be deter- 
mined by the kind of the modified grav itational theory which we are considering (for a detailed study of reconstruction 
for various modified gravities, see [l9|, |20|). The important point is that the realistic history of the expansion of the 
universe can be realized from the modified gravity reconstructed. 



III. FINITE-TIME FUTURE SINGULARITIES IN F{R) GRAVITY 

In this section, we investigate i^(i?)-gravity models and show that some models generate several known types of 
the finite-time future singularities. This phenomenon is quite natural as modified gravity represented as the Einstein 
gravity with the effective ideal fiuid with the phantom or quintessence-like EoS (see the explicit transformation in 
Ref. [21|). In some cases, it is known that such ideal (quintessence or phantom) fluid induces the finite-time future 
singularity. We present several examples which predict the accelerating FRW solutions ending at the finite-time future 
singularity by using the reconstruction technique. We demonstrate that not only the Big Rip but other three types 
of the finite-time future singularities may appear. 



As the first example, we consider the case of the Big Rip singularity [23| , where H behaves as 

i/ = -^, (3.1) 

where Hq and ig ^ire positive constants and H diverges at t = ig- To find the F{R) gravity generating the Big Rip 
singularity, we use the method of the reconstruction, namely, we construct F{R) model realizing ony given cosmology 
using the technique proposed in Ref. [20] (for the related study of reconstruction in F{R) gravity, see [23]). The action 
of F{R) gravity with general matter is given as follows: 

S^ f d^X^^{F{R) + Anattcr} • (3.2) 

The action ()3.2p can be rewritten by using proper functions P{4') and Q{(t)) of a scalar field (/) [20]: 

S= f d^Xy/^ {Pi(^)R + Q(0) + ^matter} . (3.3) 

One may regard the scalar field </> as an auxiliary scalar field because (f> has no kinetic term. By the variation over (f>, 
we obtain 

= P'{<P)R + Q'{<P) , (3.4) 

which could be solved with respect to (j) as (j> = (t>{R)- By substituting (f) = (f>{R) into the action p.3p . we obtain the 
action of F{R) gravity given by 

F{R) = P{c^{R))R + Q{(t){R)) . (3.5) 

By assuming p, p could be given by the corresponding sum of matter with a constant EoS parameters Wi and 
writing the scale factor a{t) as a = aoe^^*% where ao is a constant, one gets the second rank differential equation: 

•^ - 2^^ - 25'(^)^ + W'{4>)P{4>) + E (1 + -0 p.oao"'^'^""^c-3(i+-)^(^) . (3.6) 

If one can solve Eq. (|3.6p with respect to P{4>), the form of Q{(j)) could be found as [20] 

Q(0) = -6 {g'{cl>)f Pi<p) - 6ff'(0)^ + E p.oao '^^+"'-^e-3(i+-.).W) . (3.7) 

i 

Thus, it follows that any given history of the expansion of the universe can be realized by some specific F{R) gravity. 
Specific models unifying the sequence: the early-time acceleration, radiation/matter-dominated stage and dark energy 
epoch are constructed in Refs. [4[. Isl. l20j . 

In case of (j3.ip . if we neglect the contribution from the matter, the general solution of (j3.6p is given by 



P(<^) = F+ {to - <^r- + P (to - <t>r- , a± ^ '^° + ^ ^ ^^° ' '°^ ^^^ , (3.8) 

when /lo > 5 + 2\/6 or /iq < 5 — 2\/6 and 

P(0) ^ (to -0)-"'«+^'/^ (icos ((to - 0)ln '^° + ^^"^"'^ ) +Bsin ((to - 0) In '^° + ^^°^" ' ^ )) , (3.9) 
when 5 + 2^6 > ho > 5 — 2\/6. Using (|3.4p . (|3.5p . and (|3.7p . we find that the form of F{R) when R is large is given 

by 

F{R) ex i?i-"-/2 , (3.10) 

for ho > 5 + 2\/6 or /lo < 5 — 2\/6 case and 

F{R) ex i?(''o+i)/4 X (oscillating parts) , (3.11) 



for 5 + 2a/6 > ho>5- 2^6 case. 

Let us investigate more general singularity [24 



-/3 



H^ho{ta-t)-P . (3.12) 



Here, ho and f3 are constants, ho is assumed to be positive and t < tg because it should be for the expanding universe. 
Even for non-integer /3 < 0, some derivative of H and therefore the curvature becomes singular. We should also note 
that Ea. p.l2p tells that the scale factor a {H ~ a/ a) behaves as 

a~e 1--^ +■■■ . (3.13) 

Here, • • • expresses the regular terms. From ()3.13|) . we find that if f3 could not be any integer, the value of a, and 
therefore the value of the metric tensor, would become complex number and include the imaginary part when i > to, 
which is unphysical. This could tell that the universe could end at i = io even if /3 could be negative or less than —1. 
Since the case /3 = 1 corresponds to the Big Rip singularity, which has been investigated, we assume /3 7^ 1. 
Furthermore, since f3 ~ corresponds to dc Sitter space, which has no singularity, we assume /3 7^ 0. When /3 > 1, 
the scalar curvature R behaves as 

R - 12ij2 ^ I2hl (to ~ ty^^ . (3.14) 

On the other hand, when /3 < 1, the scalar curvature R behaves as 

R^6H ^6hof3{to-ty^~^ . (3.15) 

We may get the asymptotic solution for P when (j) -^ t^. 

• /3 > 1 case: We find the following asymptotic expression of P{4>): 

P{4,) ~ o<''"W-'»<'"-«""'(i„-*)'"'(ico8(u((„-«-'"')+B8m(a,(t„-*)-''+')) , 

When (j) -^ to, P{4>) tends to vanish. Using (|3.4p . (|3.5p and (|3.7p . we find that (at large R) F(R) is given by 

FiR) (X e^''°/'^''"''KT^) i?-i/4 X (oscillating part) . (3.17) 

• 1 > /? > case: We find the following asymptotic expression of P{4>): 

Pi^) ^ Be-("°/2(i-«)(*«-^'"' (to - cj^f^'^/' . (3.18) 

Hence F{R) is given by 

F{R) - e-iha/2(l-t3))i-6phoR)^''-^^^^^+^^j^7/8 _ ^^_^g'^ 

• /3 < case: The asymptotic expression of P{(f>) is as follows: 

P(0) ~ Ae-(''«/2(i-/3))(*o-0)^-'' (i^ _ ^yiP^-~6P+l)/8 (3_20) 

Thus F{R) is given by 

FiR) ^ (-6/io/3i?)('''+'^"'')/''''-*-'' g-(/.o/2(i-/J))(-6ho/3ii)('^-)/(''+i) _ (321) 

Note that -6/io/3i? > when ho, R> 0. 

We found the behavior of the scalar curvature R from that of H in (|3.12p . namely, when f3 > 1, R behaves as in 
p.l4p . and when /3 < 1, R behaves as in p.l5p . Now, conversely we consider the behavior of H from that of R. When 
R behaves as 

R <^ 6H ^ Ro (to - ty , (3.22) 



if 7 > 2, which corresponds to /3 = 7/2 > 1, H behaves as 



H-]J^{to-tr^\ (3.23) 



if 2 > 7 > 1, which corresponds tol>/? = 7 — 1>0, iJis given by 



^-^(7^(^o-^n" ' (3.24) 

and if 7 < 1, which corresponds to /3 = 7 — 1 < 0, one obtains 

H^Ho+ ^,^\A to-t)-^^' . (3.25) 

6(7-1) 

Here, Hq is an arbitrary constant, and it does not affect the behavior of R. Hq is chosen to vanish in (|3.12p . ^From 
H = a{t)/a{t), if 7 > 2, we find 

a{t) oc exp ((^ - 1) \[^ito tr'^'A , (3.26) 



when 2 > 7 > 1, a{t) behaves as 



a{t) (X exp ( g^^^°_^^ (to - ty ] , (3.27) 



a{t) oc exp ( Hot + — -— ^ — - {to - t) '' ] . (3.28) 



and if 7 < 1 , 

67 (7 - 1) 

In any case, there appears a sudden future singularity [25| at t = to . 

Since the second term in p.25p is smaller than the first one, one may solve p.6p asymptotically as follows: 

P'-Po(l + ^(to-0)'-') ' (^-^^^ 



with a constant Po, which gives 



F{R) - FoR + Fii?2/3/(/3+i) . (3.30) 



Here, Fq and Fi are constants. When > /3 > —1, we find 2/3/ (/3 + 1) < 0. On the other hand, when (3 < —1, we 
obtain 2/3/ (/3 + 1) > 2. As we saw in p.lOp . for f3 < —1, H diverges when t -^ to. Since we reconstruct F{R) so that 
the behavior of H could be recovered, the F{R) generates the Big Rip singularity when R is large. Thus, even if R 
is small, the F(R) generates a singularity where higher derivatives of H diverge. 

Even for F{R) gravity, we may define the effective energy density pcs and pressure p^s as (|2.3p . We now assume 
H behaves as p.l2|) . For f3 > I, when t ^ to, a ~ exp(/io (to — t) ~^ / {(3 — 1)) -^ 00 and poff, boff | ^00. If /3 = 1, 
we find a ~ (to — t) ° — > 00 and poff, \Pes\ —^ 00. If < /? < 1, a goes to a constant but p, \p\ — > 00. If —1 < /3 < 0, 
a and p vanish but |peff| -^ 00. When /3 < 0, instead of (|3.12p . as in (|3.24p . one may assume 

H^Ho + ho{to-t)-^ . (3.31) 

Hence, if — 1 < /3 < 0, poff has a finite value 3Ho/k^ in the limit t ^ to. If /3 < —1 but (3 is not any integer, a is 
finite and pcff and PeS vanish if Ho = or pcft and Poff are finite if Hq 7^ but higher derivatives of H diverge. We 
should note that the leading behavior of the scalar curvature R does not depend on Ho in (|3.3ip . and that the second 
term in (|3.3ip is relevant to the leading behavior of R. We should note, however, that Ho is relevant to the leading 
behavior of the effective energy density peff and the scale factor a. 

In Ref. J2fl|. the classification of the finite-time future singularities has been suggested in the following way: 

• Type I ("Big Rip") : For t — > t^.. a — > 00. p -^ 00 and \p\ -^ 00. This also includes the case of p, p being finite 
at Cg . 



• 



Type II ("sudden") : For t -^ tg, a -^ ag, p -^ ps and \p\ — > oo 

Type III : For t^tg^a^ras, p^oo and \p\ — > oo 

Type IV : For t -^ t^, a ^ a^, p ^ 0, |p| ^ and higher derivatives of H diverge. This also includes the case 
in which p (p) or both of p and p tend to some finite values, while higher derivatives of H diverge. 



Here, tg, 0^(7^ 0) and ps are constants. We now identify tg with io- The Type I corresponds to /3 > 1 or /3 = 1 
case. Type II to —1 < /3 < case. Type III to < /3 < 1 case, and Type IV to /3 < — 1 but /3 is not any integer 
number. When the phantom dark energy was studied, it was found that the Big Rip type (type I) singularity could 
occur. After that, it was pointed out that there could be another kind of singularity corresponding to type II [25j . 
Then, it was found in Ref. [2g that there are other kinds of singularities corresponding to type III and IV. Note 
that if only higher derivatives of Hubble rate diverge, then some combination of curvature invariants also diverges 
what leads to singularity. Thus, we have constructed several examples of F{R) gravity showing any type of the above 
finite-time future singularities. This is natural because it is known that modified gravity may lead to the effective 
phantom/quintessence phase [3| , while the phantom/quintessence-dominated universe may end up with the finite-time 
future singularity. 

The reconstruction method also tells that there appears Type I singularity for F{R) = R + aR"^ with n > 2 and 
Type HI singularity for F{R) = R — /3i?~" with n > 0, where a and /3 are constants. Note, however, that even if 
some specific model contains the finite-time future singularity, one can always make the reconstruction of the model 
in the remote past in such a way that the finite-time future singularity disappears. Usually, positive powers of the 
curvature (polynomial structure) help to make the effective quintessence/phantom phase to become transient and to 
avoid the finite-time future singularities. The corresponding examples are presented in Refs. [24. i27l|. 

IV. ABSENCE OF SINGULARITY IN MODIFIED GRAVITY 

In this section, we study various scenarios to resolve the finite-time future singularities. The most natural scenario 
is based on additional modification of the inhomogcneous EoS or the gravitational action by the term which is not 
relevant currently. In fact, however, such a term which may be relevant at very early or very late universe may resolve 
the singularity. We note that the presence of such a next-to-leading order term does not conflict with the known local 
tests. Another scenario is related with the account of quantum effects becoming relevant near to singularity. 

Let us consider the conditions for G in (j2.4p . which prevents the finite-time future singularities appearing. For 
simplicity, we only consider the case in which the matter has a constant EoS parameter w. A key equation is 
Eq. (|2.2ip . If Eq. (|2.2ip becomes inconsistent for any singularities, the singularity could not be realized. First, we 
put G = in (|2.2ip . In this case, a solution satisfying (|2.21[) is given by 



li w < —1, i.e., the phantom phase, there appears a Big Rip singularity. We may now assume u; > — 1. Even if w 
is not constant, in case —00 < w < —1, there occurs the Big Rip type singularity, where H diverges in finite future. 
For example, we may consider the case in which w behaves as w = —1 — Wq {Iq ~ t)^^ when t ^ t^. Here, Wq and 
77 are constants. In order that w could be less than — 1, wq should be positive. Hence the Hubble rate behaves as 
H ^ (to — t)~ ~^ . Thus, if 7] > —1, H diverges and the Big Rip type singularity could occur. We should note that 
when ?7 > 0, w ^ — 1 in the limit of i ^ tg. 

If H evolves as p.3ip . the r.h.s. in (|2.2ip behaves as 



_iKi^(,„_,)-2/5 when/3>I 

-^ (2iJ + 3(1 + w)H') ^ <! _ 2/3/»o+3a+^)/»g (^^ _ ^y2 o</3<l • (4.2) 

-'^ {to - t)-'^-^ 0>/3>-l 



If /3 > —1, /3 7^ 0, which correspond to Type I, II and III singularities, the l.h.s. on (|4.2p . i.e., G (H,H,- ■ ■ ] m 

Eq. (j2.2ip . diverges. One way to prevent such a singularity appearing could be that G is bounded. In this case, 
Eq. (|2.2ip becomes inconsistent with the behavior of the r.h.s. of (|4.2p . An example is as follows: 



G^G„i±||,, (43) 
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where a and b arc positive constants and Gq is a constant. 

We should also note that when /3 > 0, which corresponds to Type I or III singularity, the r.h.s. of (|4.2p becomes 
negative. Hence, if G is positive for large H, the singularity could not be realized. 

Another possibility is that if G contains the term like y/l — a'^H'^, which becomes imaginary for large H, (|4.2p could 
become inconsistent. Thus, singularities where the curvature blows up (Type I, II, III) could be prevented appearing. 
This mechanism could be applied even if w < —1. One can add an extra term as 



Gi(i7) = Go(^/l-|^-lj , (4.4) 

to G. Here, Gq and Hq are constants. If we choose Hq to be large enough, Gi is not relevant for the small curvature 
but relevant for large scale and prohibits the curvature singularity appearing. Moreover, it is easy to check the EoS 
of the effective ideal fluid induced by modified gravity. Some indications of the presence/ absence of the singularities 
can be found here. Nevertheless, one should not forget that the description of the effective ideal fluid corresponds to 
the consideration of the theory in another (Einstein) frame. The transformation of the point where the singularity 
appears from one frame to another may be tricky in modified gravity [28[. 

As the universe approaches to the singularity, the curvature may become large again. Hence, it is necessary to take 
into account quantum effects (or even those of modified gravity through effective action approach [29[). One may 
include the massless quantum effects by taking account of the contribution of the conformal anomaly as the back 
reaction near the singularity. The conformal anomaly Ta has the following well-known form [26| : 

TA = b(Fw + ^aR\+b'g + b"DR. (4.5) 

Here, Fyj is the square of 4d Weyl tensor, and it is given by 






(4.6) 



In general, with N scalar, A^i/ospinor, A^i vector fields, N2 (= or 1) gravitons and iVno higher derivative conformal 
scalars, h and h' are given by [2^ 



b 



N + &N1/2 + l2Ni + QIIN2 - 8iVHD 



120(47r) 



2 



, _ N + llAfi/2 + 62JVi + UIIN2 - 28A^HD 

360(4^)2 ■ ^^-^^ 

As is seen, & > and 6' < for the usual matter except the higher derivative conformal scalars. Notice that b" can 
be shifted by the finite renormalization of the local counter term B? , so that b" can be an arbitrary coefficient. For 
the FRW universe, we find 



Fw = , g = 24 [Hli' + i?4 j . (4.8) 

In terms of the corresponding energy density pA and pressure pA, Ta is given by [26| 

Ta = -PA + 3pA , (4.9) 

and the energy density could be conserved 

PA + 3i7(pA+PA) = 0. (4.10) 

Hence we find 

PA = -a-^ J dta^HTA , pA ^ ^ - ^ J dta^HTA . (4.11) 



By using the above expressions, Eq. (|2.3p could be modified as 



Peff + /0A = 4^'' P.a+PA^-\{2H + iH^] . (4.12) 
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If we redefine the effective energy density and pressure as 

Poff = Pcff + PA , PcS=Pctf+PA, (4-13) 

Eq. (|2.10p could be modified as 

Peff = wpcs + G (h, H, H, ■ 
G = G (f, H, • • • ) + y- - ( i + w") a-^ j dta^HTA . (4.14) 

It has been indicated [30] that the exphcit account of such quantum effects tends to moderate the singularity, to make 
the space-time less singular or at least to delay the Rip time. Clearly, these scenarios may be applied to any specific 
modified gravity under consideration. Furthermore, there remains the possibility that the action of modified gravity 
is changed by the term which is not relevant now and does not influence the local tests in such a theory currently. 
This kind of term, however, may be relevant at the very early universe or at the very late universe in such a way that 
the universe avoid s ap proaching to the singularity (in other words, say, the phantom phase becomes transient as in 
the model of Ref. ^). 

As the explicit realization of the above modification, one can consider the effective ideal fluid with the complicated 
EoS depending on the Hubble rate: 

PcS ^w(h,H,---j pcS , (4.15) 

namely, the EoS parameter could be a function oi H, H, H, ■ ■ ■ . The specific form of the dependence on the Hubble 
rate is defined by the equivalent model of modified gravity. Using (|2.3[) , we find 



= 2H + 3H^ + 3w(H,H,---]H'^ . (4.16) 

If we take iJ to be a constant Hq, Eq. (|4.16p is reduced to 

= 1 + w(Ho,H = 0,--- ) , (4.17) 



which is an algebraic equation and hence has a solution, e.g., there could be a de Sitter space solution. For example, 
if w is given by 



{H,H,...)^^ + f{H), (4.18) 



it follows 



As a special case, if w is given by 



H = hoVl-fiO) ■ (4.19) 



w[h,h)=-1 + ^, (4.20) 

Eq. (j4.16|l is trivially satisfied for any H, and therefore any cosmology can be a solution. This theory, however, has 
no predictive power. One non-trivial example is as follows: 

^i/) = -1 + 4^i/-i+V^ . (4.21) 

The solution of (|4.16p is given by the exact form of p.l2p . i.e., H = ho {to — t) . Hence, any type of the finite-time 
future singularities can be realized by the form of w in (|4.20p . Another non-trivial example is expressed by the 
following form: 
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where to J hi and hi are eonstants satisfying hi ^ hi > 0. An exact solution of (|4.16|) is given by 

i/ = ^^i±^-^^i^tanhl. (4.23) 

2 2 ^Q 

In the Umit t —f — oo, we find H -^ hi. On the other hand, in the hmit t —f c», H —^ hi. Hence the dc Sitter universe 
could be realized in both limits t -^ ±oo. Thus, we may identify the limit of i — > — oo with inflation and the limit of 
t —^ oo with the late time acceleration. Adding the effective ideal fluid with such a EoS to the model admitting the 
singularity, one can always check if this addition resolves it. 

V. FINITE-TIME FUTURE SINGULARITIES IN SCALAR-GAUSS-BONNET AND MODIFIED 

GAUSS-BONNET GRAVITIES 

In this section, we discuss the scalar-Gauss-Bonnet and modified Gauss-Bonnet gravities which predict the late-time 
acceleration ending at a finite-time future singularity. We present the corresponding effective potentials by using the 
reconstruction technique. 

Let us consider the reconstruction [19l . l3ll . |32| of the string-inspired scalar-Gauss-Bonnet gravity proposed as dark 
energy in Refs. |33l. |34|| for the investigation of the finite-time future singularities: 



j d^x^—g A_ia^0aM<^„y(0)_^(0)g 



(5.1) 



The theory depends on two scalar potentials. The explicit example motivated by string considerations is the 
following |33j : 

1/(0) = K)e-'^/*° , ^(0) = eoe'*/^« , (5.2) 

where Vq, (J)q and ^o a-re constant parameters. By choosing the parameters properly, dark energy universe ending at 
the Big Rip singularity emerges. 

There exists another example corresponding to the following choice of the potentials [34I : 

y(0) ^ \ (g^ + il^-H^X _ ^c-2^/^" - 3C/o L + ^-^c~^/A esi^/*oe9otoo^/*o ^ 

K \ to J K Iq \ to J 

Here f/Oj 9l^ and C/q are constants. The Hubble rate is given by 

i/ = 5o + ^, (/) = 0olnf. (5.4) 

Hence, when t is small, the second term in the expression of H in (|5.4p dominates and the scale factor behaves as 
a '^ t^^ . Thus, if gi — 2/3, a matter-dominated period, in which a scalar may be identified with matter, could be 
realized. On the other hand, when t is large, the first term of H in (|5.4p dominates and the Hubble rate H becomes 
constant. Consequently, the universe is asymptotically de Sitter space, which is an accelerating universe. 
In the case without matter, the equations of motion in the FRW metric are given by 

= ^H^ + -<^2 + y + 24g3 ^^^; '^ , 5.5 

K^ 2 at 

= 4 (277 + m^\ + U- - y(0) - m-'^^^ - xmn'M^ - i.n^m^ , (5.6) 

K^ V / 2 dt'^ at at 

^ 4> + 3H^ + V'{^)+^'{(j))g . (5.7) 

One may rewrite the above equations to the following form: 

hh dt uaj ox 

dt \ a dt J 



^i/ + 0^-8a^(-l-:^) . (5.8) 
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The potentials are expressed as 



K^ 2 J a[ti) \K 



^'i ^\H{t,) + ^{t,f 



ant)) 



\h''H[uY 



a(ti) /■*! dta 



^■H{t2)+Mt2f 



Equations (j5.9p show that if we consider the theory including two functions g{t) and /((/>) as 

1 



(5.9) 






2/'(0)2 "33'(/(0))e^(^(^»t/(0), 



m ^ IJ 



1 f\^ /'(-^i^'^'^^'^^ru ^ 

2— '^(<Pl) ' 



5'(/(0i))^ 



C/(0) 



rf<^i/'(0i)e- 



-9(/(0i)) _„" 



r5"(/(0i)) 



/'(0i)^ 
the solution of the field equations is given by 

</) = /-i(i) (t^fm, a = aoe^W (i? = ,g'(i)) . 



(5.10) 



(5.11) 



It is easy to include matter with constant w = Wm- In this case, it is enough to consider the theory where U{4') in 
V{(t)) and £^{(f>) in (|5.10p is replaced by 

t/(0) = J d0i/'(0i)c-^(/(^^)) l^-^g" ificf,,)) + j^ + (1 + u;„)goe-3(i+»'"WW^))) . (5.12) 

We obtain the previous solution with ag = (go/po) "'" • 

As a cousin of the scalar-Gauss-Bonnet theory, we may consider the modified Gauss-Bonnet theory [3, ISa, ISq . 
whose action is given by 



j4 / — 
a xy'—g 



R 

2k 



—2+nQ) 



(5.13) 



By introducing an auxiliary scalar field </>, we can rewrite the action (|5.13p to a form similar to the action of the scalar 
Gauss-Bonnet theory in (|5.ip : 



S = I d x^—g 



R 

2k 



^-v{^)-acf>)g 



In fact, by solving the equation for 



o = v'i<p)+c{^)g , 



(5.14) 



(5.15) 



with respect to as (/> = 4>{G) and substituting the expression of into the action (|5.14[) . we re-obtain the action (|5.13|) 
where F{G) is given by 



F{g)^-v{^{g)) + a'l^mG ■ 



(5.16) 



An example [Tg| is given by (|2.19|) . When Fq > 3/2k'^, there are two solutions, which describe effective phantom 
universe and admit the Big Rip singularity. When Fq < 3/2«;^, we have a solution, which describes the effective 
quintessence, and the other solution, which does effective phantom. 

Another example with dust, whose energy density behaves as p = poda~^, is as follows |32| : 

2C' 



VU) = ^^coth2(C0)-3CC/ocoth(C0)sinh2/^(C(A) 
^(0) == ^ / d0isinh-^/3(C0)cosh2(C0) . 



(5.17) 
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Here, C and Uq are constants. The explicit solution is given by 

a{t) = aoe^^*) , g{t) = - In (sinh (Ct)) , pod = ^^ ■ (5.18) 

Eq. (|5.18p indicates that, when (p = t is small, 5(0) behaves as g{cf>) ~ {2/3)\n(f> and, therefore, the Hubble rate 
behaves as H{t) = g'{t) ~ (2/3)/t, which surely reproduces the matter-dominated phase. On the other hand, when 
(j) ^ t is large, g{(j)) behaves as g ^ (2/3) ((70), namely, H ^^ 2C/3 and the universe asymptotically goes to de Sitter 
space. Hence, the model given by (j5.17p with matter shows the transition from the matter-dominated phase to the 
accelerating universe, which is asymptotically de Sitter space. 

We may consider the reconstruction of F{Q) gravity [l^[34|. The field equations in the FRW background are given 

by 

-- ^H2^y(0) + 24ij3^iM!)) , (5.19) 

at 



K 



Wozv , o^2\ ,.,.^ o^^d'a^t)) ^.r^i^d^m)) ,.r.^dam) 



K 



= -^{2H + 3H^)- V{(t)) - 8H^ ^y^ " - 16HH ^'y " - 16H'' ^y^ " , (5.20) 



df^ dt dt 



which could be rewritten as 



im))^\fdt,^^w{ty, vm))^^H(tf-Mt)H{t)w{t), w{t)^^f^^H{ty. (5.21) 

Since there is no kinetic term of 0, there is a freedom to redefine <j) as (p -^ ip ~ 'f{4')- By using the redefinition, we 
may choose the scalar field as a time coordinate: (j) = t. Thus one gets 

!/(</>) = -g'(0)'- 35' (0)e«(*)t/(0), m^^J ^<^i^(^^('^i) ' ^('^)^^y d0ie-^(^^V'(0i) • (5.22) 

The solution is given by 

a = aoe9(*) {H = g'{t)) . (5.23) 

As in the case of the scalar-Gauss-Bonnct theory, we may include matter. 

Let us apply the reconstruction program to the study of the finite-time future singularities, where the Hubble rate 
behaves as 

H^hoit.-ty" . (5.24) 

First, we consider the modified Gauss-Bonnet gravity. For (3—1 and ha > case, namely, usual Big Rip singularity, 
since g'{(j)) = H{(j)), we find 

5(0) = -/ioln^V^, (5.25) 

to 

where to is a constant of the integration. C/(0) in (|5.22p has the following form: 



Uicf,)^^ I rf0W:i— ^li)'" ^" ., =1^" ,i^#^^'^ '^^'"'^ -TT' : ■ (5-26) 



2 r., r ts - 01 y" ho I Uo~ ^.tAfL-i) (i- " '^) ^^^^ ^0 ^ 1 

K- .1 \ to J {ts-M^ l^ln^ when/io = l 



Here, C/q and ti are constants of the integration. At the next step, one gets 

3''oto"'^^o , 3hl{h„ + l) , , , , 



^ 3^(l-21n|) when/io = l 



^0 - 8h-i3-ho) + a.-{ho~i) ^hen ho y 1,3 
a^) - { -^Inf + tC^ when/io = 3 • (5.28) 

^o-f^lnf + 3)$' when/^o = f 
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Here, ^0 Eind ^2 arc integration constants but they are irrelevant to the action because the Gauss-Bonnet invariant is 
a total derivative and ^o and ^2 correspond to the constant shift of the coefficient of the Gauss-Bonnet invariant. We 
further redefine the scalar field (^ as $ = t, — (p. The obtained forms of F($) and ^($) do not contain the parameter 
ts corresponding to the Big Rip time. Hence, the resulting F(Q) does not contain tg either and tg could be determined 
dynamically by initial conditions. 
When f3 ^ 0, g{(f>) is given by 



gW 



/3-1 



[ts 



.1-/3 



-50 



(5.29) 



where go is a constant of the integration but the constant is irrelevant and does not appear in the final expressions of 
V{(j)) and ^{(j)). We therefore choose go = 0. The Uicf)) has the following form: 



C/(0) 



2hoP 



hiU-M-'-'e-^^^^-^^^'~' 



2ho(3 



«M/3-l) 



dxa;i/(^-i)e-/^^ 



where 



X = {ts- (f)) 



1-/3 



(5.30) 



(5.31) 



We now consider the case /3 > 1, which corresponds to the Type I singularity. In this case, a; — > (X) when — > i^ 
When X is large, the following expression can be used: 



X 

1 — 5-a; 

a a^ 



a ^„_i a{a-l) 2 , 



1 



(3-1 



Keeping the leading order, one finds 



C/(0) 



2/3 



-e f^ 



and therefore 



3h. 



K^{ts 



6/3^-/3-1 3/^2 



1^($) ii$-^'^ + -^<P-''-' ^^-'f , ^($) 



'^jits 



-2/3 



$2/3 
8^ 



-Co 



Here, ^ = tg — (jj again and ^o is an irrelevant constant of the integration. 

In case of /3 < 1, which corresponds to the Type II, HI, IV singularities, we find that a; ^ when 
the expression 



1 



r.a+1 



„a+2 



r.a+3 



a + 1 



a + 2 



2! (a + 3) 



and only keeping the leading order, we find 



t/(0) = ^ {t. 



Uo, 



(5.32) 

(5.33) 

(5.34) 

ts. Using 

(5.35) 

(5.36) 



where Uq is a constant of the integration. Hence, 
y($) ^ _£f|>ci>-2/3 _ 3/,q[/„$-/3 ^ ^(ci,) . 



$ 



0+1 



t^olf^W 



4K^ho{/3+l) 8h;^(2/3+l) 
-1 

1 — 



1 In ± ^ '^0^ 



2K^ho sTil t2 



^0 when /3 7^ -1, -i 
when /3 = — 1 
when /3 = — -i 



(5.37) 



Here, ^0 and ^2 are (irrelevant) constants of the integration. 

In case of the scalar Gauss-Bonnet gravity, there is a freedom or ambiguity in the choice of f{4>)- For simplicity, 
we now choose 



t = /(0) = «:2^ 



(5.38) 
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The equations 


(|5.inD 


arc rewritten as 












V{v) 


- i^'(^)^ 


1 
2k4 








e(^) 


= Ifd^, 


qsM 
9'iViY 








u{v) 


= / dipic^ 


V' 


where ip ^ k'^(J) 


. By the calculation similar to the case 


oiF{g 


Rip), 


we find 











-3g'M)e«('^)t/(^) 
-^9 (<<5i) + - 



(5.39) 



C/(0) 



f^o 



Iho 



riy (i^ - "P) ' 



In 



K2to7r5(?l(,_l) 

ts-v (t.-y)''"+^ 






when ho y^ 1 
when /iQ = 1 



Here, Uq and ti are constants of the integration again. Subsequently, 

2^(^+1) ^hen /lo 



F($) = 



^hO + l 

3 



2$2 



1 



21nf 

to 



1 



t';;°uo'i>'^-''° 



^($) 



^ to" Co* " , $2 , $* U J. / 1 Q 

?0 - 8h^(3-fto) + 8K2(h.o-l) + 32h-'a(ho + l)K* ^^^^ /iQ 7^ i, 3 



72 ^^h 



Co 



iln^ 

2 ^^^ ti 



16k2 



$2_L ** 



32h2(^o + l)K* 



when Hq — 3 
when ft-o = 1 



Here, $ = ig — iy9. In terms of $, the kinetic term of is rewritten as 



1 






In case of /3 > 1 , one obtains 

y($) 
In case of /3 < 1, we get 

C/(^) 
F($) 



K 



W __^(t,_^)i-'3 



'^) 



-e z^- 



(t3-V3)''_^_ JiO 



/iqk" 



3 



:g^<i>-2^-^, ?($) 



$2/3 

3^ 



-eo- 



$3/3+1 

8/i3 (3;9 + 1) k4 



2 (ts-^) '' + t/0- 



2^0 

3_^2 



i:!!0$-2/3 
2 



* 



/a+i 



C('i>) 



"4k2/,o(/3+1) 

1 In* 






^$1-/3 

k4 



4k2/jo 

' 2n'^ho 



S,hj{2fj+l) 



t2 ' 8/i^ 

^In^ 



-I- 

,23 + 2 



16/tj5K'4(/3+l) 



16/i2k''(/3+1) 



*^ 



le/il^K-iC/a+i) 



when /3 7^ — 1, 
when j3 = —1 
when /3 = — i 



(5.40) 



(5.41) 



(5.42) 



(5.43) 



(5.44) 



(5.45) 



Thus, any type of the finite-time future singularities can be realized in F(Q) gravity and the scalar-Gauss-Bonnet 
one. We have constructed the specific examples of the above models containing the finite-time future singularity 
by using the reconstruction program. We should note, however, that the relation between the appearance of the 
singularity and the forms of V{4>) and ^{4>) in the scalar-Gauss-Bonnet gravity is not so clear. In other words, even if 
the explicit forms of T^((/>) and ^(^) are given, it is difficult to understand if the theory could generate the singularity 
or not. The corresponding investigation of the asymptotic behavior of solution sh ould be d one in order to answer to 
this question. In case of F{Q) gravity, however, if F{Q) contains the term like \/Gq — G^ (Go > 0), which becomes 
imaginary, and therefore inconsistent, if \G\ > Go-: the singularity should not appear. In other words, even if the 
solution contains the finite-time future singularity, the additional modification of the action may resolve it. 
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VI. NON-MINIMAL MAXWELL-EINSTEIN GRAVITY 

111 this section, we consider some cosmological effects in the non-minimal Maxwell-Einstein gravity with general 
gravitational coupling. We account for our model and derive the effective energy density and pressure of the universe. 

A. Model 

We consider the following model action [33 : 

•S'gr = d'^xy/^[CE}i + Cem] , (6.1) 

Ceh = T^R, (6.2) 

Cem = -\l{R)Ff,,F^''' , (6.3) 

I{R) = 1 + I{R) , (6.4) 

where F^,^ = 9^/1^ — 9^/1^ is the electromagnetic field-strength tensor. Here, A^ is the U{1) gauge field. Furthermore, 
I{R) is an arbitrary function of R. It is known that the coupling between the scalar curvature and the Lagrangian 
of the electroma gnet ic field arises in curved space-time due to one-loop vacuum-polarization effects in Quantum 
Electrodynamics [3^. (In Ref. [39|, a non-minimal gravitational Yang-Mills (YM) theory, in which the YM field 
couples to a function of the scalar curvature, has been discussed. Note that this study maybe generalized also for 
F{R) gravity coupled to non-linear electrodynamics [ioj.') 

Taking variations of the action Eq. (|6.ip with respect to the metric 5^^ and the U{1) gauge field A^, we obtain the 
gravitational field equation and the equation of motion of Af^ as 37 



R^. - Ig^^.R - «;'r(E^^ , (6.5) 



with 



4' 



2 

and 



^ll'{R)Fc.pF^I'R^, + g^^U [l' {R)F^pF''f] - V^V. [l' {R)F^pF^I'] \ , (6.6) 



^ a^(^/^/(i?)i^^=0, (6.7) 



respectively, where Tyi^ is the contribution to the energy-momentum tensor from the electromagnetic field. 

B. Effective energy density and pressure of the universe 

We now assume the fiat FRW space-time with the metric in Eq. (|2.1[) . We here consider the case in which there 
exist only magnetic fields and hence electric fields are negligible. In addition, only one component of B is non-zero 
and hence other two components are zero. In this case, it follows from divB = that the off-diagonal components of 
the last term on the r.h.s. of Eq. ()6.6|) for T^j, , i.e., V^Vi/ [f'{R)FapF°'^^^ are zero. Thus, all of the off-diagonal 

components of T^^, are zero (for the argument about the problem of off-diagonal components of electromagnetic 
energy-momentum tensor in non- minimal Maxwell-gravity theory, see [4l|). Moreover, because we assume that there 
exist the magnetic fields as background quantities at the 0th order, the magnetic fields do not have the dependence 
on the space components x. 

In the FRW background, the equation of motion for the U{\) gauge field in the Coulomb gauge, d^Aj{t,x) ~ 0, 
and the case of Ao(t, x) = 0, becomes 

• 1 (3) 

A,{t, x)+{H+^\ A(^t, x)-—AA,{t,x)^0, (6.8) 
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(3) 

where A = d^di is the flat 3-dimensional Laplacian. It follows from Eq. (|6.8p that the Fourier mode Ai{k,t) satisfies 
the equation 



A,{k,t) +iH+-\ A,{k,t) + ^A,{k,t) = 0. 

Replacing the independent variable t by conformal time rj = j dt/a{t), we find that Eq. (|6.9p becomes 

d^A,{k,rj) , 1 dI{7^)dA{k,Tj) 



drf 



I{jl) drj drj 



+ k^Ai{k,ri) = 0. 



(6.9) 



(6.10) 



It is impossible to obtain the exact solution of Eq. (|6.10p for the generic evolution of the coupling function / at the 
inflationary stage. However, by using the WKB approximation on subhorizon scales and the long-wavelength approx- 
imation on superhorizon scales, and matching these solutions at the horizon crossing [42| . we find an approximate 
solution as 



\A{k,v)\' = \C{k)\' 



1 



2kl{i]k) 



1- 



1 1 rf/(%) 

2 kl{i]k) d-q 



+ i 



k r '-^dh 



(6.11) 



where r]k and rjf are the conformal time at the horizon-crossing and one at the end of inflation, respectively. From 
Eq. (|6.1ip . we obtain the amplitude of the proper magnetic fields in the position space 



IB. 



(proper) 



it)? 



fc|C(fc)|2fc4 



(6.12) 



on a comoving scale L = 27r/fc. Thus, from Eq. (|6.12p we see that the proper magnetic fields evolves as |iJj(P™P'='^' (t)|2 = 
|i?P/a^, where \B\ is a constant. (The validity of this behavior of the proper magnetic fields, namely, that \B\ is a 
constant, is shown in Appendix.) This means that the influence of the coupling function / on the value of the proper 
magnetic fields exists only during inflation. (On the other hand, because the expression of the energy density of 
the magnetic fields is described by that of the magnetic fields multiplying I due to the Lagrangian (|6.3p , the energy 
density of the magnetic fields depends on / also after inflation. We can see this point from the first term on the r.h.s. 
of Eq. (j6.13p shown below.) The conductivity of the universe CTc is negligibly small during inflation, because there are 
few charged particles at that time. After the reheating stage, a number of charged particles are produced, so that 
the conductivity immediately jumps to a large value: CTc ^ H. Consequently, for a large enough conductivity at the 
reheating stage, the proper magnetic fields evolve in proportion to a~^{t) in the radiation-dominated stage and the 
subsequent matter-dominated stage |43| . 

In this case, it follows from Eq. (j6.6p that the quantity corresponding to the effective energy density of the universe 
Poff and that corresponding to the effective pressure poff arc given by 



Pes 

PcS 



HR) , 



- ( 5i/2 + h) I'{R) + QH [AHH + H) I"{R) 



\B\ 



-H^ + 5H) r(R) - 6 ( -2QH^H + 4H- - HH + H] I" [R) 



\B\ 



zqUhh + h] I"'{R) 



(6.13) 



(6.14) 



where we have used the following relations by taking into account that electric fields are negligible: g'^'^FopFoa — 
(l/4)gooi^a/3F"/' = |i?,(P™P^')(i)|V2, and F^^^F'^'' = 2|i?/P™P^')(i)|2. 

Finally, we remark the following point. Suppose that /(i?) is (almost) constant at the present time. We now assume 
that for the small curvature, I(R) behaves as 



1{R) ~ /oi?" 



(6.15) 



with constants /q and a. Here, we consider the case a < 0. The energy density of the magnetic fields is given by 
PB = (1/2) \B.,^^^°^''^\t)\^I(R) = [|B|V (2a4)] I(R). Here we take de Sitter background as the future universe. In 
such a case, when R tends to zero in the future, the energy density of the magnetic field becomes larger and larger 
in comparison with its current value. Hence, if for the small curvature, non-minimal gravitational coupling of the 
electromagnetic fields behaves as ~ /q^" with a < in the future, the strength of current magnetic fields of the 
universe may evolve to very large values in the future universe. 
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VII. FINITE-TIME FUTURE SINGULARITIES IN NON-MINIMAL MAXWELL-EINSTEIN GRAVITY 

In this section, we investigate the finite-time future singularities in non-minimal Maxwell-Einstein gravity. We 
study the forms of the non-minimal gravitational coupling of the electromagnetic field generating the finite-time 
future singularities and the general conditions for the non-minimal gravitational coupling of the electromagnetic field 
in order that the finite-time future singularities cannot emerge. 

It follows from (|2.2p that the FRW equations are given by 

\h^ ^ Pes, (7.1) 

--^(2iJ + 3i/2) = Pes, (7.2) 

where poff and Pcff sue given by Eqs. (j6.13p and (|6.14p . respectively. 

We investigate the form of I{R) which produces the Big Rip singularity, 

i/~-^, (7.3) 

to - 1 

where ho is a positive constant, and H diverges at t ~ tQ. In this case, 

12/ig + 6ho .^ 

H ^ — "2— , a -- ao (to - Ej , (7-4) 

(to — t) 

where ao is a constant. We now assume that for the large curvature, /(-R) behaves as Eq. (|6.15p . Hence peff in 
Eq. (|6.13p behaves as (to — t)" ""*" ", but the l.h.s. in the FRW equation iH"^ / k^ = poff does as (to — t)~ . The 
consistency gives 



namely 



Eq. (|7.ip also shows 



-2a + 4/io, (7.5) 



a — 1 

ha = ——— or a = 1 + 2/io . (7.6) 



^^l = h {12KI + 6hoy-' I (^^''o+J^'^o) + 3 ^_^ ^^2hl + 6ho) {ho + bKl) +&a{a- 1) ho {2ho + 4/.g)] I M 



Av2 




(7.7) 

which requires that /g should be negative. In the second line of (|7.7p . we have used (|7.6p and deleted a. 

As a result, it follows from Eqs. (|6.15p and (|7.6I1 that the Big Rip singularity in Eq. (j7.3p can emerge only when 
for the large curvature, I{R) behaves as i?^+^''°. If the form of I{R) is given by other terms, the Big Rip singularity 
cannot emerge. We here note that if exactly I(R) = loR", H = ho/ (to — t) is an exact solution. 

Next, we study the form of I{R) which gives a more general singularity in Eq. (|3.12p . In this case. 



R^Qho 



/3 + 2ho (to - t)-(^-^)l {to - t)-(^+^) , a ^ ao exp I ^ [to t)"'^^"^' 



(7.8) 



We also assume that for the large curvature, I{R) behaves as Eq. (|6.15p . li (3 < —1, in the limit t ^ to, R ^ 0. 
Hence, we consider this case later. If /3 > 1, a —> oo, and hence pcs -^ and pcs -^ because pcs oc a~^ and 
Pcs oc a~*. On the other hand, H -^ oo. Thus Eqs. (|7.ip and (|7.2p cannot be satisfied. 

If a > and < /3 < 1, pes in Eq. ((6T3l) evolves as (to - t)""^''+^\ but the l.h.s. of Eq. dLj) does as (to - t)"^*^. 
Thus, the consistency gives 

-2/3 = -a(/3+l) , (7.9) 



namely, 



From Eq. (|7.1|) . we also find 
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a 2(3 , , 

(3=- or a^-^. 7.10 

2 — a p + I 



Zhl _ /q (6/10/?)" (1 - P) l^p 
«2 2a4(/? + l) ' ^'-''^ 

where we have used Eq. (|7.10p . and on the l.h.s. we have taken only the leading term. Eq. (|7.1ip requires that /q 
should be negative. Consequently, if a > and < /? < 1, in the limit t ^^ t^, a ^ a^, R ^t oo^ Poff -^ oo, and 
Ipcff I — * oo. Hence the Type III singularity emerges. If a > and — 1 < /3 < 0, p^s -^ oo, but H ^ 0. Hence Eq. (|7.1|) 
cannot be satisfied. 

If (/9 — 1) / (/3 + 1) < a < and — 1 < /3 < 0, in the limit t ^ to, a ^ Qq, R — > oo, Poff -^ 0, and \pcs\ -^ oo- 
Although the final value of poff is not a finite one but vanishes, this singularity can be considered to the Type II. The 
reason is as follows. In this case, when / and H are given by / = 1 + IqR" and H = Hq + h^ {to — t)~ , where Hq is a 
constant, respectively, in the above limit poff — > Po- From Eqs. (|6.13p and (|7.1|1 . we find po = SHq/k^ = \B\^ / {2a^. 
Hence, po is a finite value. 

If a < (/3 — 1) / (/3 + 1) and —1 < /3 < 0, in the limit t —t t^, a ^ ao, R —>■ oo, pcff ^ 0, and |pcff | -^ 0, but i7 — > oo. 
Hence Eq. (|7.2p cannot be satisfied. If a < and < /? < 1, pcff ^ 0, but H ^ oo. Hence Eq. (|7.ip cannot be 
satisfied. 

In addition, we investigate the case in which /3 < — 1. In this case, in the limit t —^ t^, a —^ uq and i? — > 0. We 
assume that for the small curvature, R behaves as Eq. (|6.15p . If a > (/? — 1) / (/? + 1), in the limit t — > to, pcff -^ 0, 
boff I — * 0, and higher derivatives of H diverge. Hence the Type IV singularity emerges. If < a < (/? — 1) / (/? + 1), 
PcS — >■ and \pcff\ -^ oo. However, H ^ and H ^ Q. Hence Eq. (|7.2p cannot be satisfied. 

We note that if I{R) is a constant (the case in which I{R) = 1 corresponds to the ordinary Maxwell theory), any 
singularity cannot emerge. 

Here we mention the case in which I{R) is given by the Hu-Sawicki form [14| 

C2 {R/m^) + 1 

which satisfies the conditions: Ivchr-^oo Ius{R) = ci/c2 = const and liTciR-tG Ius{R) = 0. Here, ci and C2 are dimen- 
sionless constants, n is a positive constant, and m denotes a mass scale. The following form [5| also has the same 
features: 

I{R)^I^o{R) = ^T . ,. —TTTT^ (7-13) 



C3 +C4 



{[{R/M^) - {RJM^)f''+^ + (i?e/M2)2«+i| ' 



which satisfies the conditions: limfl_>oo ^no (^) = l/c4 = const and limfl^o^No(^) = 0. Here, C3 and C4 are dimen- 
sionless constants, g is a positive integer, M denotes a mass scale, and Re is current curvature. If /3 < — 1 and I{R) 
is given by /hs(^) in Eq. (|7.12p or /no(^) in Eq- (|7.13p . in the limit t ^ to, a^ ao, R^ 0, pcS -^ 0, and \pcb\ -^ 0. 
In addition, higher derivatives of H diverge. Thus the Type IV singularity emerges. 

Consequently, it is demonstrated that Maxwell theory which is coupled non-minimally with Einstein gravity may 
produce finite-time singularities in future, depending on the form of non-minimal gravitational coupling. 

The general conditions for I{R) in order that the finite-time future singularities whose form is given by Eqs. (|7.3p or 
(|3.12p cannot emerge are that in the limit t -^ to, I{R) -^ I, where I{y^ 0) is a finite constant, I'{R) -^ 0, I"{R) -^ 0, 
and r"{R) -> 0. 

VIII. INFLUENCE OF NON-MINIMAL GRAVITATIONAL COUPLING ON THE FINITE-TIME 
FUTURE SINGULARITIES IN MODIFIED GRAVITY 

In this section, we consider the case in which there exist the finite-time future singularities in modified F(R) 
gravity and investigate the influence of non-minimal gravitational coupling on them. We show that a non-minimal 
gravitational coupling of the electromagnetic field can remove the finite-time future singularities or make the singularity 
stronger (or weaker). 
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In this case, the total energy density and pressure of the universe are given by ptot = PcS+Pmg andptot = PcS+PmG: 
respectively. Here, Peff and p^g are given by Eqs. (|6.13p and (|6.14p . respectively. Moreover, it follows from Eqs. 
and (|2.9p that pmg and puG are given by 



1 

PMG = — 
hi 



Pmg — — 

hi 



1 



f{R) + 3(h^ + h) f(R) - 18 [AH^H + HH] ,f"{R) 



\f{R) - i/iH'^ + h) f'{R) + 6 Uh'^H + AH^ + QHH + 'h\ f"{R) + 36 UhH + h\ f"'{R) 



(8.1) 
(8.2) 



In this case, it follows from Eqs. ((2?2)) . (|6A3)) . ((6ll)) . ((81|) . and (|8?2)) that the FRW equations arc given by 



—^H — Ptot — 



— {2H + 3H' 



_(m 



\~ 



-(5H^ + h) I'{R) + QH [AHH + H) I"{R) 



\B\ 



Ptot 



- {F{R) -R) + 3(h^ + h) {F'{R) - 1) - 18 [AH^H + HH] F"{R) 



liR) 



(8.3) 



H^ + 5H) I'{R) - 6 -20H^H + AH^ - HH + H I"{R) 



-36 [ah H + H^ I"'{R) 



\B[ 



^-{F{R)-R)-[3H^ + H]{F'{R)-l) 



+6{8H^H + AH'^ + 6HH + h) F"{R) + 36{AHH + H) F"'{R) 



(8.4) 



Using Eqs. ([13]) and dH]), we find 
-I{R) 







2(-8H^ + H]l'{R) + 6 (32H'^H - AH'^ + AHH - H] I"{R) - 36 (AHH + H] I"'{R) 



\B[ 



2HF'{R) + 6 ( -AH^H + AH^ + 3HH + H ) F"{R) + 36 UHH + H] F"'{R) 



(8.5) 



In modified gravity with the ordinary Maxwell theory, if F(i?) is given by Eq. p.30p . i.e., F(R) ~ FgR+FiR"^, where 
g is a constant, a finite-time future singularity appears. Let us account for non- minimal gravitational electromagnetic 
theory in this F(i?)-gravity model. 

When H behaves as 



H^ho (to - t)" 



(8.6) 



where u is an positive integer, there exists no finite-time future singularity. In what follows, we consider the case 
u>2. From Eq. dM]), we find 



R^6H 



-6uho (to — t) 



«-i 



a ^ ao exp 



^° (to-t)"+^' 



u+1 



■7) 



where in the expression of R we have taken only the leading term. 

We investigate the form of the non-minimal gravitational coupling of the electromagnetic field I{R) which produces 
the solution (|8.7p . We here assume that I{R) behaves as (j6.15p . The first and second terms on the r.h.s. of 
Eq. (j8.5p are the non-minimal gravitational electromagnetic coupling and the modified-gravity sectors, respectively. 
When t is close to to, the leading term of the non-minimal gravitational electromagnetic coupling sector evolves 
as (to — ty^~ '^°'~ ' . On the other hand, if g < 1, or (7 > 1 and u < q/ {q — 2), that of the modified-gravity 



sector does as (to — t) 



(«-l)(g-l)-2 



Hence the consistency gives a ~ q. Moreover, in order that the leading term 



of the non- minimal gravitational electromagnetic coupling sector should not diverge in the limit t ^ to , a must be 
a > (u + 1) / (u — 1). Taking only the leading terms in Eq. (|8.5p and using a ~ q, wc find 



lo 



4Fi 

|B|2k2 



(8.8) 



If g > 1 and u > (7/(5 — 2), the leading term of the modified-gravity sector behaves as (to — t)^ 



consistency gives a = 2u/ (m — 1). In this case, taking only the leading terms in Eq. 



Hence the 
and using a — 2u/ {u — 1), 
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we obtain 



/o 



a^^o 



\B\^K^6u^iu + l) 



-6hou) 



-2/(u-l) 



(8.9) 



Consequently, we see that the non-minimal gravitational coupling of the electromagnetic field I{R) ^ /o-R" with the 
specific values of Iq and a stated above can resolve the finite-time future singularities which occur in pure modified 
gravity. 

Next, we study the case in which non-minimal gravitational coupling of the electromagnetic field does not remove 
the singularity but makes it stronger (or weaker). We also assume that for the large curvature, I{R) behaves as (|6.15p . 
Using the result in Eq. (j3.10|) . we consider the case in which for the large curvature, F{R) behaves as F{R) oc R^, 
where q = 1 — a^/2 < 1. In this case, the Big Rip singularity in Eq. (|7.3p emerges. It follows from Eq. (|7.6p that 
if a = 1 + 2/io, in the limit t ^ to, PcS ^ oo and |pcff | ~^ oo. Hence, the non-minimal gravitational coupling of the 
electromagnetic field makes the singularity stronger. 

When there exist a more general singularity in Eq. (|3.12p with < /3 < 1, which is the Type III singularity and 
can appear for the form of F{R) in Eq. (j3.19p . and a > 0, in the limit t — > to, Pcff ^ oo arid \pcs\ -^ oo. Thus the 
non-minimal gravitational coupling of the electromagnetic field makes the singularity stronger. If —1 < /3 < 0, namely, 
there exists the Type II singularity, which can appear for the form oi F{R) in Eq. (|3.2ip . and (/3 — 1) / (/? 4- 1) < a < 0, 
Poff —* and |poff| -^ oo. For |ptot| > |pmg| (|ptot| < |pmg|), therefore, the non-minimal gravitational coupling of 
the electromagnetic field makes the singularity stronger (weaker). Thus, the non-minimal gravitational coupling in 
Maxwell theory may qualitatively infiuence to the universe future. For instances, for some forms of non-minimal 
gravitational coupling it resolves the finite-time future singularity or it may change its properties. 



IX. CONCLUSION 

In the present paper, we have considered the finite-time future singularities in modified gravity: F{R) gravity, 
scalar-Gauss-Bonnct or modified Gauss-Bonnet one and effective fluid with the inhomogcneous EoS. It has been 
demonstrated that depending on the specific form of the model under consideration, the universe may approach to 
the finite-time future singularity of all four known types. It is not easy to say from the very beginning if the particular 
theory brings the accelerating universe into the singularity or not. As a rule, each theory should be checked to see if 
the singularity appears on the FRW accelerating solutions. It is also interesting that the finite-time future singularity 
cannot be seen through the study of local tests: theory passing known local and cosmological tests may produce the 
future universe with (or without) the singularity. We have presented several examples of modified gravity having 
the accelerating dark energy solutions with the finite-time future singularity by explicitly using the reconstruction 
program. Our results are shown in the Table U below, where the possibility of four types of the finite-time future 
singularity is explained for modified gravities under consideration. 



Moreover, we have discussed several theoretical scenarios resolving the finite-time future singularity. In particular, 
we have considered the additional modification of the fluid with the inhomogcneous EoS or that of the gravity action 



TABLE I: Summary of the results: po, c, ci, C2, C3, ci, £2, and C3 are proper constants. 





Type I 

(/3>1) 


Type II 
(-l</3<0) 


Type III 
(0 < /3 < 1) 


Type IV 
(/3 < -1, /3 : not integer) 


Einstein gravity with 
inhomogeneous EoS 


p + pocp('3+l)/2/5 


/ 1/2 l/2\(/'+l)//3 
P + poc (pi/2 ^ p^/ j 


F{R) gravity 


F{R)oce'^^~^R-i 




F(il) oc e-'^-R^^ R3 


F(i?)oce-=-«''+'R 8(/3+i) 


Scalar-Gauss-Bonnet 
gravity 


y(*) oc^-^", 5($) (x$2'' 


C ci<I.'^+-^+C2-I'^'^+^ (/3^ -1,-1/2) 
[ ci<I>i/2 + c2ln'I> {13 = -1/2) 


Modified Gauss-Bonnet 
gravity 


Vi<S>) oc*-=^, 


y(*) oc ci^-'^l^ + C2^-f + C3*^-'^, 
( ci^./'+l + C2*2/3 + l + c3*2/3+2 (^ ^ _1_ „l/2) 
^(<I>) OC .^ ciln* + C2*-^-|-C3*2^+2 (/3 = -l) 

\ ci*l/2 + 02 In* -1- C3<I>2''+2 (/3 = -1/2) 


Non-minimal Maxwell- 
Einstein gravity 


I{R) oc R'-+'-"^o (/3 = 1), 
Not possible if /3 > 1. 


Possible 


Not possible 


Possible 
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by the term which is relevant at the very early (or very late) universe. If the corresponding term is negligible at 
the current epoch, such modification is always possible, at least, from the theoretical point of view. Nevertheless, in 
order to check if the corresponding term is realistic, its role at the very early/late universe should be confirmed by 
using observational data. Another scenario to remove the singularity is related with the account of quantum effects 
(or even quantum gravity). However, the corresponding consideration may give the preliminary results, at best, due 
to the absence of the consistent theory of quantum gravity. 

Furthermore, we have investigated the non-minimal Maxwell-Einstein (or Maxwell-_F(i?)) gravity in the similar 
fashion. Wc have studied the forms of the non-minimal gravitational coupling generating the finite-time future sin- 
gularities and the general conditions for the non-minimal gravitational coupling in order that the finite-time future 
singularities cannot emerge. In addition, we have considered the influence of the non-minimal gravitational coupling 
on the finite-time future singularities in modified gravity. As a result, it has been shown that the non-minimal gravi- 
tational coupling in Maxwell-modified gravity can remove the finite-time future singularities or make the singularity 
stronger (or weaker). 

The observational data of the current dark energy epoch still cannot tell what is its exact nature: phantom, ACDM 
or quintessence type. It has been demonstrated in this paper that in some modified gravities with the effective 
phantom or quintessence EoS the finite-time future singularity can emerge. In this respect, the interpretation of 
the observational data confirming (or excluding) the approach to the finite-time future singularity is fundamentally 
important. From one side, it may clarify the distant future of our universe. From another side, it may help to define 
the evolution of the universe and the current value of the effective EoS parameter: how close is w to —1? 
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APPENDIX: 

We here consider Eq. (|6.7p . Let us now assume 

E, = Fo^ = doA - d,Ao = , (A.l) 

where Ei is the electric field. Then 

A, = 9, /" dtAo + a . (A.2) 

Here Ci does not depend on time t. Then we find 

Fy = d,Aj - djA, = d,C, - djC, . (A.3) 

Hence, the Fij and therefore the magnetic flux Bi = eijkFjk docs not depend on time. Since the metric g^i^ and hence 
the scalar curvature R do not depend on the spatial coordinates, Eq. (|6.7p reduces to the following form: 

d'F.j = ACj - djd'a = . (A.4) 

The solution of (|A.4p is given by the constant Fij or constant magnetic flux Bi. 
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